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ABSTRACT 

In this paper, we study the photon-photon pair production optical depth 
for gamma-rays with energies from 10 to 200 GeV emitted by powerful blazars 
due to the diffuse radiation field of broad line region (BLR). There are four key 
parameters in the BLR model employed to determine the 7—7 attenuation optical 
depth of these gamma-rays. They are the gamma-ray emitting radius FL,, the 
BLR luminosity Lblr, the BLR half thickness h and the ratio tblr//cov of the 
Thomson optical depth to the covering factor of BLR. For FSRQs, on average, it 
is impossible for gamma-rays with energies from 10 to 200 GeV to escape from 
the diffuse radiation field of the BLR. If GLAST could detect these gamma- 
rays for most of FSRQs, the gamma-ray emitting region is likely to be outside 
the cavity formed by the BLR. Otherwise, the emitting region is likely to be 
inside the BLR cavity. As examples, we estimate the photon-photon absorption 
optical depth of gamma-rays with energies from 10 to 200 GeV for two powerful 
blazars, HFSRQ PKS 0405-123 and FSRQ 3C 279. Comparing our results with 
GLAST observations in the future could test whether the model employed and 
the relevant assumptions in this paper are reliable and reasonable, and then limit 
constraints on the position of the gamma-ray emitting region relative to the BLR 
and the properties of the BLR. 

Subject headings: galaxies: active — galaxies: jets — gamma rays: theory - 
quasars: individual (3C 279, PKS 0405-123) 



1 National Astronomical Observatories/ Yunnan Astronomical Observatory, the Chinese Academy of Sci- 
ences, Kunming, Yunnan 650011, P. R. China. 

2 Graduate School of the Chinese Academy of Sciences, Beijing, P. R. China; send offprint requests to 
liuhongtaollll@hotmail.com 



- 2 - 



1. INTRODUCTION 

Thanks to the development of gamma-ray astronomy, the sample of GeV or TeV gamma- 
ray emitters is increasing. With the detection of high energy gamma-rays in over 60 blazars, 
such as 3C 279 and PKS 1510-089, in the GeV energy range by the EGRET experiment 
aboard the Compton Gamma Ray Observatory (Catanese et al. 1997; Fichtel et al. 1994; 
Lin et al. 1997; Mukherjee et al. 1997; Thompson et al. 1995, 1996; Villata et al. 1997), an 
exceptional opportunity presents itself for the understanding of the central engine operating 
in blazars. Some of blazars have been also firmly detected by atmospheric Cherenkov tele- 
scopes (ACTs) at energies above 1 TeV, such as Mrk 421 (Punch et al. 1992) and Mrk 501 
(Quinn et al. 1996). The Large Area Telescope instrument on the Gamma- Ray Large Area 
Space Telescope (GLAST), the next-generation high energy gamma-ray telescope which will 
be launched next year, will observe gamma-rays with energies from 20 MeV to greater than 
300 GeV, and have the unique capability to detect thousands of gamma-ray blazars to red- 
shifts of at least z — 4, with sufficient angular resolution to allow identification of a large 
fraction of their optical counterparts (see e.g. Chen, Reyes, & Ritz 2004). The High Energy 
Stereoscopic System (HESS) is an imaging atmospheric Cherenkov detector dedicated to 
the ground based observation of gamma-rays at energies above 100 GeV (Funk et al. 2004; 
Hinton 2004; Hofmann 2003), and successfully detected TeV gamma-rays from many sources 
(see e.g. Aharonian et al. 2005). The Major Atmospheric Gamma Imaging Cherenkov 
telescope (MAGIC) is, currently the largest single-dish Imaging Air Cherenkov Telescope in 
operation, the one designed to have the lowest energy threshold of ~ 30 GeV among the 
new Cherenkov telescopes (see e.g. Baixeras et al. 2004). The MAGIC telescope detected 
TeV gamma-rays from HESS J1813+178 (Albert et al. 2006). GLAST, combined with 
new-generation TeV instruments such as VERITAS, HESS and MAGIC, will tremendously 
improve blazar spectra studies, filling in the band from 20 MeV to 10 TeV with high signif- 
icance data for hundreds of active galactic nuclei (AGNs) (see Gehrels & Michelson 1999). 
Future measurements of the gamma-ray spectra shape and its variability of blazars would 
tremendously improve our understanding to blazars. 

Blazars are radio-loud AGNs characterized by luminous nonthermal continuum emission 
extending from radio up to GeV/TeV energies, very rapid variability, high and variable 
polarization, flat radio spectrum, high brightness temperatures and supperluminal motion 
of compact radio cores (e.g. Urry & Padovani 1995). These unusual characteristics of blazars 
are generally believed to be attributed by emission from a relativistic jet oriented at a small 
angle to the line of sight (Blandford & Rees 1978). The spectral energy distributions of 
blazars consist of two broad bumps. The first component is from the synchrotron radiation 
processes. The second component is contributed by inverse Compton emission of the same 
electron population, including synchrotron self-Compton (SSC) scattering synchrotron seed 
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photons, and external Compton (EC) scattering seed photons from sources outside the jet 
(see e.g. Bottcher 1999). The external soft photon fields not only provide target photons for 
EC processes to produce the second component, but also absorb these gamma-rays from the 
EC processes, because gamma-rays with energies above 10 GeV interact with optical-UV 
photons and then could be attenuated by photon-photon pair production collisions with the 
external soft photons. The positions of gamma-ray emitting regions are still an open and 
controversial issue in the researches on blazars. It is suggested that gamma-rays are produced 
within broad line region (BLR) and that the gamma-ray emitting radius i? 7 ranges roughly 
between 0.03 and 0.3 pc (Ghisellini & Madau 1996). It is argued by Georganopoulos et al. 
(2001) that the radiative plasma in relativistic jets of powerful blazars are within cavities 
formed by BLRs. However, it is also argued by other researchers that the gamma-ray emitting 
regions are outside BLRs (Lindfors, Valtaoja, & Turler 2005; Sokolov & Marscher 2005). 

Observations found a large population of X-ray strong FSRQs, labeled HFSRQs by 
Perlman et al. (1998). HFSRQs have strong broad emission lines, flat hard X-ray spec- 
tra similar to classical FSRQs, but high synchrotron peak frequencies and steep soft X-ray 
spectra similar to intermediate or high synchrotron peak frequency BL Lac objects (e.g., 
Georganopoulos 2000; Padovani et al. 2003). Therefore, it is expected that HFSRQs are 
likely to be GeV-TeV gamma-ray emitters. For example, RGB J1629+4008 has a syn- 
chrotron peak ~ 2 x 10 16 Hz and a predicted Compton flux vF v of nearly 10~ 10 erg s _1 cm -2 
at 10 25 Hz (Padovani et al. 2002). However, the BLRs can reprocess and scatter the central 
UV radiation to produce diffuse line emission and continuum radiation, respectively. The 
two external diffuse radiation could absorb gamma-rays with energies from 10 to 200 GeV. 
Then, it is unknown whether these gamma-rays could be detected by GLAST even if blazars 
intrinsically produce gamma-rays with energies from 10 to 200 GeV in the observer's frame. 
In this paper, we focus our efforts on the photon-photon absorption for these gamma-rays 
by the two diffuse radiation fields of the BLRs in the powerful blazars, and then hope to 
limit constraints on the positions of gamma-ray emitting regions by comparing our results 
with GLAST observations in the near future. 

The structure of this paper is as follows. Section 2 consists of four subsections, in which 
subsection 2.1 presents the cross section of pair production, subsection 2.2 the photon- 
photon attenuation optical depth, subsection 2.3 the diffuse radiation field of the BLR, and 
subsection 2.4 the dependence of r 77 on parameters. Section 3 presents the photon-photon 
absorption optical depth for FSRQs, section 4 presents applications to individual sources, 
and the last section is for discussions and conclusions. Throughout this paper, we use a flat 
cosmology with a deceleration factor q = 0.5 and a Hubble constant H = 75 km s -1 Mpc -1 . 
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2. THEORETICAL CALCULATION 

2.1. The Cross Section of Pair Production 

Gamma-ray photons can be absorbed by the photo-annihilation/pair creation process 
7 + 7 — >• e + + e~, for which the total cross section is (e.g., Gould & Schreder 1967; Jauch & 
Rohrlich 1976): 



(3_^ ln i±£ _2x(2-x 2 ) 
1 — x 



(1) 



where <tt is the Thomson cross section, and a; is the velocity of the pairs in the center of 
moment system and defined as 



where h is the Plank constant, m e the rest mass of electron, c the speed of light, e± (e 2 ) is 
the energy of the high (low) energy photon in units of the rest mass energy of electron, and 
/j, is the cosine of the collision angle 6 (Coppi & Blandford 1990). The threshold condition 
of photon-photon pair creation process is 

jUie 2 (l -cos0) > 1. (3) 

The total cross section reaches the maximum value with cr max ~ 0.26o"t at x ~ 0.7, which 
corresponds to 

£l£2(l - COS0) R3 4. (4) 

Then we have £i£ 2 ^ 2 because < 6 < it. For the 7 — 7 attenuation process, only the 
photons with the product e\S2 > 2 could have the maximum cross section <j max ~ 0.26o" T . 



2.2. The 7 — 7 Attenuation Optical Depth 

The 7 — 7 attenuation optical depth per unit length dl at distance R from the central 
black hole along the direction of jet axis is 



d_ 
Jl 



r jy (e u R) = J J J <T(x)n ph (R,e 2 ,ti)(l - n)de 2 dQ, 



(5) 



where n p h(R, e 2 , f2) is the number density of soft photon at the position R per unit frequency 
range per unit solid angle. For these gamma-rays produced at the position of i? 7 along the 
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jet axis, the 7 — 7 attenuation optical depth is 

f°° d 

r 77 (£i) = J — t 77 (£i, R)dR, (6) 

and the escape probability of gamma photons of £ 1; -P e sc( £ i), is 

P csc (e 1 ) = e~ T ^\ (7) 

In next paragraphs, we calculate the energy density distribution of gamma-ray absorbing 
photons from both scattering and reprocessing, and then estimate the 7 — 7 attenuation 
optical depth. 

First, we calculate the intensity by integrating the emissivity along lines of sight through 
BLR, which is assumed as a spherical shell of clouds with the outer radius rBLR,out and the 
inner radius rBLR.in (see Fig. 1). A fraction df cov (r) = n c (r)a c (r)dr of the central UV 
luminosity Luv would then be reprocessed into emission lines in r — > r + dr, where n c (r) 
and cr c (r) are the number density and the cross section of clouds at radius r, respectively. 
Since a thin spherical shell has volume 4irr 2 dr, the emissivity (in erg s _1 cm -3 ster" 1 ) of 
reprocessed radiation inside the BLR at the radius r is 

/ x _ -^uv^/cov _ L\jvn c cr c 

Jhney ) — 77 097 — ~~rp, 9 9 5 V->) 

where a c = nR 2 and -R c is the radius of clouds at the radius r. For n c = ^ c o( r / r BLR,in)~ p and 
R c = -R c o( r / r BLR,m)' J , where n c o and R c q is the number density and the radius of clouds at 
?"BLR,m, respectively, the covering factor f cov of the entire BLR to the central UV radiation 
is 

/•r-BLR.out / j. \ 2 1~P 

fcov = / n c0 irR 2 c0 dr. (9) 

•/r-BLR.in V^BLR.in/ 

In this paper, we adopt the power law exponents preferred by Kaspi & Netzer (1999), q — 1/3 
and p = 1.5. The central UV radiation is also scattered by the BLR clouds into the diffuse 
continuum, and the diffuse soft photons can absorb gamma-rays with energies above 10 
GeV. A fraction g?tblr — n c aj:dr of the central UV luminosity would then be scattered in 
r — > r + dr, where n e is the number density of electrons in the BLR clouds at the radius r. 
Since this spherical shell has volume 4iir 2 dr, the emissivity of scattered radiation inside the 
BLR at the radius r is 

/ x _ -^uv^tblr _ L V yn e a T , 1f ,x 

if the Thomson optical depth of the BLR tblr *C 1. For n e = n e o(^/TBLR,m) _s , where n e o is 
the number density of electrons in the BLR clouds at rBLR,in, the Thomson optical depth of 
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the entire BLR to the central UV radiation is 

/'''BLR, out / r \ S 

tblr = / «co^t dr. (11) 

JrBLR.in V^BLR.in/ 

In this paper, we adopt the power law exponent preferred by Kaspi & Netzer (1999), s — 1. 
From equations (8) and (9) and the relation of the luminosity of broad emission lines from 
the BLR and the central UV luminosity Lblr = /cov-^uv (D'Elia et al. 2003), we get the 
expression for the diffuse line emission 

' ( \ = -^blr^ 2<? ~ p ~ 2 / 19 x 

3]ine{f) - ^ 2 rr-BLR,out r . 2q -p dr - 

-''"BLR, in 

From equations (10) and (11) and Lblr = /cov^uv, we also get the emissivity for the diffuse 
continuum radiation 

/ \ _ L BLR T BLR T~ S ~ /.ox 

The intensity of radiation with angle 9 to the jet axis at a position R is given by 

^max 

I(R,9)= / j(r)dl, (14) 
Jo 

where j(r) represents j cont (r), jn ne (r) or j cont (r) +j Vmc (r), and vanishes when r > r BL R,out or 
r < r BLR ^ in . We can get the following relations from the triangle in Figure 1 

r 2 = R 2 -2Rlcos6 + l 2 , (15) 

l max = Rcos6 + rJ cos 2 6 + y - i. (16) 

The energy density of radiation from the BLR is expressed as 

U(R) = J J ll^p-dtt=' 2 ^ jT I(R,9)sm9d9, (17) 

and is a function of distance R from the center along the direction of the jet axis. 



2.3. The Diffuse Radiation Field of BLR 

Though, the emission lines Lya, CIV, Mgll, if 7 , Hp and H a are the major contributors 
to the total broad line emission, the contributions from all the emission lines reported by 
Francis et al. (1991) and the emission line H a reported by Gaskell et al. (1981) are included 
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to estimate the total flux of broad emission lines, F B lr- Together 35 components of emission 
lines are considered to estimate the total flux of broad emission lines. The sum of line ratios 
N u , relative to the Lya with a reference value of N„ = 100, of all the 35 components of 
emission lines is equal to 555.77 (Francis et al. 1991). The relative intensities of the 35 
components of emission lines are presented in Figure 3. It is assumed that the flux of a 
certain emission line is equal to 



BLR- 



li 



555.77 

According to the definition of flux, we have relations of Fblr = / I\ine(R, &) cos 8d£l and 
F v = J I(R, u, 8) cos 8dQ. Comparing the two relations with equation (18), we have a re- 
lation of I(R,u,8) = N u I\ ine (R, #)/555.77 between the total intensity Iu ne (R,8) and the 
monochromatic intensity I(R, v, 8) for the emission line. Then the photon number density 
of the emission line along the direction of 8 at the position of R is 

nph(i? ' Q) = = 555.77 hue ■ (19) 

A diluted blackbody spectrum at the position of R along the direction of 8 is assumed for 
the diffuse continuum as 

n ph (R, v, O) = gr^y 71 = < R i #Kb(^, T )> (2°) 

where k is the Boltzmann constant, and n(R, 8) is the normalization factor to a blackbody 
spectrum ri\,\,{u,T) and can be calculated by the formula 

<m= J cont{R,e ] — , (2i) 

c J u l n hh (v,T)hvdv 

because I cont (R,8) = cn(R,8) J n^iv ^T^hvdv . In the calculations of the 7 — 7 attenuation 
optical depth, the soft photon frequency is considered in the range from v% = 10 14 ' 6 Hz 
to i>2 = 10 16 ' 5 Hz for the diffuse continuum at the bands of the optical to UV, and the 
blackbody temperature T = 10 5 K is assumed for the inner regions of the accretion disk. 



2.4. The Dependence of r 77 on Parameters 

According to the interpretation of the data obtained in many reverberation campaigns, 
production of broad emission lines in quasars is peaked around some distance (see Peterson 
1993; Kaspi 2000; Sulentic et al. 2000). Wu et al. (2004) derived an empirical relation 
between the BLR size tblr deduced from the reverberation mapping and Hp luminosity 

logr BL R= -1-695 + 0.684 log— ^ pc. (22) 

10 42 erg s -1 
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From equation (18) and the line ratio of N u = 22 for the emission line ifg, the empirical 
relation is also expressed as 



where Lblr,45 = -^blr/(10 45 ergs -1 ). We assume that the distance tblr deduced from 
equation (23) represents the median of r B LR,in and r BLR,out> that is, r B LR,m = r blr — h and 
^BLR,out = ^blr + h, where h is the half thickness of the spherical shell of clouds. 

We derived an average broad line luminosity of Lblr = 10 45 ' 35 erg s -1 for 25 FSRQs 
from Cao & Jiang (1999), and took it as the typical value of Lblr for powerful blazars. The 
typical value of Lblr = 10 45 35 erg s _1 for FSRQs corresponds to a distance r blr — 0.43 pc. 
In order to illustrate the radial dependence of the energy density of the diffuse radiation 
on the two parameters tblr and h, equivalently rsLR.in and rBLR,out, we firstly consider four 
cases with four values of h at a fixed r blr = 0.43 pc, and then consider four cases with four 
values of r blr at a fixed h = 0.20 pc by assuming a ratio Tblr//cov = 1- The relevant energy 
density distributions of the diffuse radiation fields are plotted in Figure 2. It can be seen that 
for each case the energy density peaks around the inner radius r B LR,in, and decreases rapidly 
beyond the outer radius ?"BLR,out- On the whole, the energy density inside the inner radius 
increases as the parameter h increases for a fixed tblr, and the average energy density inside 
the outer radius also increases as h increases. The energy density inside the inner radius 
decreases as the parameter t B lr increases for a fixed h, and the average inside the outer 
radius also decreases as v blr increases. The energy density outside the inner radius does not 
monotonously increase with h for a fixed tblr (see Fig. 2c). At a fixed tblr, the increasing 
of h results in the increasing of rBLR,out and decreasing of rBLR.in- It can be seen from 
Figure 2c that the variation of the inner radius can significantly change the energy density 
distributions, but the energy density distributions slightly alter as the outer radius changes. 
At a fixed h, the increasing of tblr results in the increasing of rBLR,in and rBLr,out- It is 
obvious that the energy density tends to become smaller as either rBLR.in or ?"BLR,out becomes 
larger (see Fig. 2a, b). The energy density outside the inner radius does not monotonously 
decrease with t B lr for a fixed h (see Fig. 2a, b). In addition, the radial dependencies of the 
energy densities have a very similar behavior as R is large enough. In order to study the 
dependence of r 77 on parameters h, i? 7 and tblr, a ratio of tblr// CO v = 1 is assumed and 
the BLR luminosity of L B lr = 10 45 35 erg s -1 is adopted for FSRQs. 

We assume i? 7 = r blr and take four values of h to study the dependence of the photon- 
photon absorption optical depth on the parameter h. The calculated results are presented 
in Figure 4. The parameter h influences the photon-photon absorption optical depth, but 
the influence of h on the calculated results is not very important (see Fig. 4a). It can 




0.684 




(23) 
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be seen from Figure 4a that the optical depth r 77 does not monotonously increase with h, 
qualitatively changing the dependence on h with E y . We study the dependence of r 77 on h 
for four values of E 1 to check if this effect is real. The calculated results are presented in 
Figure 4b, and we can see from Figure 4b that the optical depth r 77 does not monotonously 
increase with h, qualitatively changing the dependence on h with E 1 . To better show this 
effect, we normalize r 77 (/i) by the factor r 77 (0.05 pc) for the four curves in Figure 4b, and 
the normalized results are presented in Figure 4c. For the curves I, II and III with the 
relevant E 1 at 50, 100 and 150 GeV, the optical depth r 77 at four h have an order of 
r 77 (0.05 pc) < r 77 (0.30 pc) < r 77 (0.20 pc) < r 77 (0.15 pc) (see Fig. 4c). For the curve 
I+II with Ej at 200 GeV, the optical depth r 77 at four h have a sequence of r 77 (0.30 pc) < 
r 77 (0.05 pc) < r 77 (0.20 pc) < r 77 (0.15 pc) (see Fig. 4c). It is obvious that the optical depth 
r 77 does not monotonously increase with h and the dependence of r 77 on h qualitatively 
changes with E 1 (see Fig. 4c). Then this effect is real, and it could result from several 
factors. Though, the energy density inside the inner radius and the average energy density 
increase as h increases for a fixed tblr, a majority of the diffuse radiation cannot make a 
contribution to the optical depth r 77 because of the assumption of i? 7 = r blr- Only the 
diffuse radiation outside v blr can make a contribution to the optical depth r 77 , and then this 
results in the non-monotonous dependence of r 77 on h at a fixed tblr because the energy 
density distributions outside tblr do not monotonously increase with h. Only the fraction 
of soft photons matching the threshold condition can make a contribution to the photon- 
photon absorption, and this makes a contribution to the non-monotonous dependence on h. 
These factors together result in this effect that the optical depth r 77 does not monotonously 
increase with h, qualitatively changing the dependence on h with E 7 . If i? 7 = r B LR,in, the 
dependence on h at a fixed tblr could be monotonous because more absorption with h 
increasing can be contributed by the diffuse radiation between rsLR,in and tblr, and this 
is justified by the photon-photon absorption for HFSRQ PKS 0405—123 in Figure 8 (see 
section 4). 

In order to study the dependence of r 77 on the parameter i? 7 , we adopt tblr = 0.43 pc 
and assume the half thickness of BLR h = 0.20 pc. The calculated results with four values 
of FLy are presented in Figure 5. The farther the gamma-ray emitting region from the 
central black hole, the less the absorption optical depth of gamma-rays. When the gamma- 
ray emitting region is close to the outer radius of BLR, the two diffuse radiation fields 
we concerned are transparent for gamma-rays with energies from 10 to 200 GeV. If the 
gamma-ray emitting region is close to the inner radius of BLR, the two diffuse radiation 
fields we concerned are not transparent for gamma-rays between 10 and 200 GeV. This 
trend of the absorption optical depth in Figure 5 can be easily explained from the energy 
density distributions in Figure 2. The closer the gamma-ray emitting region to the center of 
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BLR, the larger the displacement of the escaping gamma-rays and the more the gamma-ray 
absorbing photons. The absorption optical depth r 77 rapidly decreases as the parameter i? 7 
increases. The photon field energy density sharply reduces as R 7 > r B LR,m, and this makes 
the absorption optical depth r 77 rapidly reduce. 

We take four values of tblr to study the dependence of the absorption optical depth r 77 
on the parameter tblr- In the calculations, it is supposed that i? 7 = tblr and h = 0.20 pc. 
The calculated results are presented in Figure 6. It can be seen in Figure 6 that the optical 
depth r 77 slowly reduces as the parameter tblr increases. This trend of the absorption 
optical depth can be also easily explained from the energy density distributions in Figure 2. 
The increasing of BLR scale tblr makes the inner and outer radii of BLR become larger, 
and this results in the decreasing of the energy density of radiation fields we concerned, and 
then the decreasing of the optical depth of photon-photon absorption. 



3. THE PHOTON-PHOTON ABSORPTION OPTICAL DEPTH FOR 

FSRQs 

According to Ghisellini et al. (1998), for FSRQs, the Lorentz factors have the typi- 
cal value of T ~ 15, the magnetic energy densities are on the order of C/ B ~ 0.1 erg cm~ 3 
measured in the blob comoving frame, and the Compton dominance of L EC /L Syn is ~ 20, 
where L EC is the EC luminosity and Ls yn the synchrotron luminosity. It is suggested that 
there are relations of Lec -^ssc and LEc/^Syn — 20 for FSRQs on average, where 
Lssc is the SSC luminosity (Georganopoulos et al. 2001). This Compton dominance 
Z/Ec/^Syn — 20 is consistent with the result of Ghisellini et al. (1998). There is a relation- 
ship of L^c/Lsyn — U' h /U B = C/ph/C/e between the ratio of EC to synchrotron luminosity 
and the ratio of external radiation energy density C/ p h to magnetic energy density Ub for 
powerful blazars. The external radiation energy densities of FSRQs have a typical value of 
C/p h ~ U B L EC /L Syn ~ r~ 2 f/ B L E c/^s y n ~ 10 -2 - 1 erg cm" 3 . Assuming the ratio of T BLR /f cov 
on the order of unity, the calculations indicate that the average of the external radiation 
energy densities within the outer radius r B LR,out is on the order of C/ ph ~ 10 -2 - 1 erg cm~ 3 as 
?"BLR,m ~ 0.2 pc and ?"BLR,out ~ 0.6 pc. The median of the probable positions of inner and 
outer radii is around 0.4 pc, which is in well agreement with tblr — 0.43 pc derived from 
equation (23) and the typical value of Lblr — 10 45 ' 35 erg s~ x for FSRQs. This agreement 
implies it reasonable that we adopt the value of tblr derived from equation (23) as the 
median of rBLR,in and ?"BLR,out for FSRQs. The typical values of Lblr = 10 45 ' 35 erg s -1 and 
C/ph ~ 10 -21 erg cm" 3 for FSRQs limit the probable positions of inner and outer radii of 
BLR around 0.2 pc and 0.6 pc, respectively, and also confine the half thickness of BLR about 



-li- 



ft, ~ 0.20 pc. Then, we adopt r B LR — 0.43 pc and h ~ 0.20 pc to estimate the photon-photon 
absorption optical depth of gamma-rays with energies from 10 to 200 GeV for FSRQs. 

However, the gamma-ray emitting regions are still a controversial issue in the studies 
on blazars. It has been suggested that gamma-rays are produced within the BLR, and the 
gamma-ray emitting radius i? 7 ranges roughly between 0.03 and 0.3 pc (Ghisellini & Madau 
1996). The radiative plasma in the relativistic jets of the powerful blazars are within the 
BLR cavities (Georganopoulos et al. 2001), and then there are R 7 < ?" B LR,m for FSRQs 
with the typical value of L B lr = 10 5 erg s _1 . Based on the dependence of r 77 on the 
parameter i? 7 discussed in previous section, the closer the gamma-ray emitting regions to 
the central black holes, the larger the optical depth of photon-photon absorption, and then 
we adopt i? 7 ~ r BLRin to estimate the 7 — 7 absorption optical depth of gamma-rays with 
energies from 10 to 200 GeV. Assuming for FSRQs that i? 7 ~ r BLRiin , h ~ 0.20 pc and 
£blr = 10 45 - 35 erg s _1 corresponding to r BL R — 0.43 pc, we estimate the photon-photon 
absorption optical depth of gamma-rays between 10 and 200 GeV. The line emission is 
completely transparent to gamma-rays below 30 GeV (see Fig. 7), because the threshold 
condition of photon-photon pair creation process cannot be satisfied between these gamma- 
rays below 30 GeV and the emission lines listed in Figure 3. However, the diluted blackbody 
radiation is not transparent to these gamma-rays below 30 GeV (see Fig. 7). On the other 
hand, the emission lines contribute much more than the diluted blackbody radiation to the 
absorption for gamma-rays with energies from ~50 to 200 GeV. Even if the diluted blackbody 
radiation is not considered in the calculations, the line emission is not transparent to these 
gamma-rays with energies from ~50 to 200 GeV (see Fig. 7). Then, gamma-rays with 
energies from 10 to 200 GeV cannot escape from the emission line and diluted blackbody 
radiation fields of the BLRs in FSRQs as long as i? 7 < ?"BLR,m (see Fig. 7). In this section, 
the redshift effect is not considered for gamma-rays with energies from 10 to 200 GeV in the 
observed spectrum, and this effect will be considered for individual sources in the following 
section. 



4. APPLICATION TO INDIVIDUAL SOURCES 

The discovery of a large population of X-ray strong FSRQs (labeled HFSRQs by Perlman 
et al. 1998) gives a fundamental change in our perception of the broadband spectra of FSRQs. 
HFSRQs have strong broad emission lines, flat hard X-ray spectra similar to classical FSRQs, 
but high synchrotron peak frequency and steep soft X-ray spectra similar to intermediate 
or high synchrotron peak frequency BL Lac objects, and are classified as a new subclass of 
blazars (e.g., Georganopoulos 2000; Padovani et al. 2003). An example of such an object 
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is RGB J1629+4008 with a synchrotron peak ~ 2x 10 16 Hz and a predicted flux vF v of 
nearly 10~ 10 erg s" 1 cm" 2 at 10 25 Hz (Padovani et al. 2002). Therefore, HFSRQs are likely 
to be GeV-TeV gamma-ray emitters with the Compton peak energies possibly from multi- 
ten to multi-hundred GeV. A powerful HFSRQ with Lblr obtained in literature is PKS 
0405—123, and it has a BLR luminosity Lblr = 10 45 ' 67 erg s" 1 and redshift z = 0.574 (e.g., 
Georganopoulos 2000). Assuming that HFSRQs have the same properties of BLRs as the 
classical FSRQs, we estimate the optical depth of photon-photon absorption for the strong- 
lined, powerful HFSRQ PKS 0405—123. In the calculations, we assume i? 7 ~ r BLRiin and 
tblr//cov = 1, and take three values of h, because it is unclear whether the parameter h 
equals ~ 0.20 pc for PKS 0405—123. From Figure 8 we see that the two diffuse radiation fields 
of BLR are not transparent to gamma-rays with energies from 10 to 200 GeV in the observed 
spectrum for this powerful HFSRQ as long as _R 7 < rBLR,in- If GLAST detect gamma-rays 
with energies from 10 to 200 GeV emitted by PKS 0405—123, either the gamma-ray emitting 
regions are not within the BLR cavity or the model used and the relevant assumptions made 
in this paper are incorrect. 

The classical FSRQ 3C 279, z = 0.536, is one of the brightest extragalactic objects in 
the gamma-ray sky, and has a BLR luminosity Lblr = 10 44 ' 41 erg s" 1 (Cao & Jiang 1999). 
It was detected by the EGRET in the 0.1-10 GeV energy domain, and its spectrum above 
~ 100 MeV can be fitted by a power law and do not show any signature of gamma-ray 
absorption by pair production up to ~ 10 GeV (Fichtel et al. 1994; von Montigny et al. 
1995). We calculate the optical depth to photon-photon absorption for 3C 279 by adopting 
r BLR,in = 0.1 pc and r BLR,out = 0.4 pc from Hartman et al. (2001) and assuming i? 7 ~ r BLR,in- 
The photon-photon absorption optical depth is presented by the dashed curve in Figure 8. 
It can be seen in Figure 8 that the diffuse radiation fields are not transparent to gamma-rays 
with energies around 10 GeV, and this result is not consistent with the observations (e.g., 
Fichtel et al. 1994; von Montigny et al. 1995; Ghisellini & Madau 1996; Wehrle et al. 1998). 
There are two probable factors that make the optical depth overestimated. The first one is 
the assumption of i? 7 ~ TBLR.m that can lead to the overestimation of r 77 . The second one 
is the assumption of tblr/ /cov = 1 that can also result in the same consequence as the first 
one. If Tblr//cov = 1/3, the photon-photon absorption optical depth of these gamma-rays 
around 10 GeV is less than unity, and the two diffuse radiation fields are optically thick to 
these gamma-rays in the interval between multi-ten and 200 GeV. GLAST observations in 
the future could limit strong constraints on the position of the gamma-ray emitting region 
relative to the inner radius of BLR and the BLR properties for 3C 279. 
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5. DISCUSSIONS AND CONCLUSIONS 

It can be seen from equation (12) that the emissivity of emission lines is unrelated to 
the covering factor, and then the relevant absorption for gamma-rays is also unrelated with 
the value of / cov . The gamma-ray absorption contributed by the diluted blackbody radiation 
is in proportion to the ratio t B lr// C ov (see eq. 13). There is a large uncertainty in the 
BLR covering factor, and the uncertainty could affect our results. Early estimates of this 
quantity indicated a covering factor f cov ~ 5 — 10%, while recent observations indicated 
instead / cov ~ 30% (e.g., Maiolino et al. 2001). Estimates of the BLR covering factors need 
to measure UV continuum emission from accretion disks. In the case of blazars, these UV 
continuum emission are swamped by the beaming components from relativistic jets, and this 
makes it difficult to calculate the BLR covering factors for blazars. The way out to roughly 
estimate f cov is to know the BLR luminosity Lblr of blazars, in which the thermal emission 
is clearly dominating the synchrotron ones. D'Elia et al. (2003) found only two of these 
blazars in literatures: 3C 273 and PKS 2149-306. They obtained f cov ~ 7% for the first 
source, while / cov ~ 10% for the second one. D'Elia et al. (2003) adopted the BLR covering 
factor f cov ~ 10% as the typical value of f cov for their blazar sample. Donea & Protheroe 
(2003) suggested f cov ~ 3% for the spherical distribution of the BLR clouds in quasars. The 
issue of the Thomson optical depth of BLRs has been rarely studied and has gotten nowhere 
for blazars. Blandford & Levison (1995) adopted the Thomson optical depth of BLR as 
Tblr — 0.01. Then the ratio Tblr//cov of the Thomson optical depth to the covering factor 
is likely to be typically of tblr//cov ~ 1/3 for the spherical distribution of the BLR clouds in 
blazars. This value of tblr/Zcov ~ 1/3 is on the same order of magnitude as tblr//cov = 1, 
assumed in the previous sections. 

As stated in section 3, the external radiation energy densities of FSRQs are on the order 
of Uph ~ 10 -2 - 1 erg cm~ 3 , and this typical value limits constraints on the probable positions 
of inner and outer radii of BLR or the half thickness of spherical shell of BLR. The probable 
positions of the inner and outer radii are around r B LR,in ~ 0.2 pc and r BLR,out ~ 0.6 pc, 
respectively. The median of this probable inner and outer radii is around 0.4 pc, which is 
consistent with the radius tblr — 0.43 pc derived from the reverberation mapping. This 
agreement implies that the radius tblr derived from the reverberation mapping probably 
denotes the median of rBLR,in and rBLR.out, and our assumption of tblr in the previous 
sections is appropriate. The energy density distribution for the BLR with rBLR,m ~ 0.2 pc, 
t BLR,out ~ 0.6 pc and r BLR / f cov ~ 1/3 has an average of energy densities within the outer 
radius on the order of U p h ~ 10~ 21 erg cm" 3 . The corresponding absorption optical depth 
for gamma-rays is presented by the curve III in Figure 7, from which it can be seen that 
the two diffuse radiation fields of the BLRs are not transparent to these gamma-rays with 
energies from 10 to 200 GeV. 
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Because of the relativistic beaming effect, most of gamma-rays for FSRQs are contained 
within a radiation cone with a half open angle of Atp ~ 1/r ~ 3.8°. If the central UV photons 
directly from accretion disks travel through the radiation cone, the UV photons can have 
photon-photon pair creation processes with gamma-rays within the radiation cone. Assuming 
the radius of the UV radiation region is around Rd ~ 100i? s ~ 0.001 pc (i? s = 2GM/c 2 
and M = 1O 8 M ), and the gamma-ray emitting region i? 7 ~ 0.1 pc, the angle between the 
jet direction and the travelling direction of UV photons at i? 7 is ~ arcsin R D /Rj ~ 0.57°. 
Then the photons within the radiation cone have the collision angles of 9 < 4.4°. For 
the UV photons at the frequency v ~ 10 16 - 5 Hz with energies of e 2 c± 2.56 x 10~ 4 and 
these gamma-rays of e 1 ~ 3.91 x 10 5 corresponding to energies around 200 GeV, the left 
of equation (3) has a upper limit of < 0.15, which is less than the value at the right of 
equation (3), and then the two kinds of photons cannot be absorbed by the photon-photon 
pair creation processes. Therefore, the central UV radiation have negligible contributions to 
the absorption for gamma-rays relative to the diffuse radiation from the BLRs. 

The energies carried by the annihilated gamma-rays would be deposited in the lower 
frequency radiation field by the synchrotron and inverse Compton processes of the electron- 
positron pairs and the pair cascade and pair annihilation processes (see Protheroe & Stanev 
1993; Sauge & Henri 2004; Zdziarski & Coppi 1991). The intense creation of pairs would 
produce a strong radiation at low energy X-rays or at GeV energy (Protheroe & Stanev 
1993; Sauge & Henri 2004; Zdziarski & Coppi 1991). The electron-positron pair cascade 
can cause the soft X-ray excesses (Zdziarski & Coppi 1991). The produced electron-positron 
pairs could make a difference at lower energies, around 1 GeV, where the pileup of photons 
below the absorption feature results in a significant flattening in the observed spectrum 
relative to the emission spectrum (Protheroe & Stanev 1993). Then the electron-positron 
pairs produced by the photon-photon annihilation processes in FSRQs could result in the 
similar observational effects mentioned above at around 1 GeV gamma-rays or low energy 
X-rays. 

In summary, we used a BLR model to study the photon-photon absorption optical 
depth for gamma-rays with energies from 10 to 200 GeV in the observed spectrum. There 
are four key parameters for the BLR model: the gamma-ray emitting radius R^, the BLR 
luminosity L BLR , the half thickness of BLR h and the ratio tblr// C ov The BLR luminosity 
Z/blr determines the distance tblr, which limits the BLR inner and outer radii of rBLR,in and 
^BLR,out if the parameter h is given. We discussed the dependence of the optical depth r 77 on 
the parameters i? 7 , h and v blr- The diffuse external radiation energy density in the gamma- 
ray emitting region is typically on the order of £7 p h ~ 10~ 2 erg cm~ 3 for FSRQs, and the 
ratio tblr/ /cov might be typically on the order of tblr/ /cov ~ 1 for blazars. Then the inner 
and outer radii of BLR are required to be around 

^BLR,in ~ 0.2 pc and TBLR,out ~ 0.6 pc, 
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respectively. It is impossible for gamma-rays from 10 to 200 GeV emitted by FSRQs to 
escape from the diffuse radiation fields of the BLRs if i? 7 ~ r B LR,in (see Fig. 7). If GLAST 
could detect gamma-rays with energies from 10 to 200 GeV for most of FSRQs, there would 
be _R 7 > r BL R,in for FSRQs on average, otherwise, R 1 < r BL R,in- If t B lr//cov ~ 1/3 and 
~ ^BLR,in for the classical FSRQ 3C 279, the diffuse soft photon field is not transparent 
to gamma-rays with energies from multi-ten to 200 GeV and is transparent to gamma- 
rays around 10 GeV. GLAST observations in the energy domain of multi-ten to 200 GeV 
could limit constraints on the position of the gamma-ray emitting region relative to the 
inner radius of BLR and the BLR properties for 3C 279. For the powerful HFSRQ PKS 
0405—123, the diffuse photon field is not transparent to gamma-rays with energies from 
10 to 200 GeV in the observed spectrum. GLAST observations in the future could limit 
constraints on the gamma-ray emitting region and the BLR in this HFSRQ. Our results are 
model dependent, especially dependent on the gamma-ray emitting radius R y and the inner 
radius rBLR.m- To determine the relationship of the inner radius rBLR,in with the radius tblr 
derived from the reverberation mapping needs more theoretical and observational researches. 
GLAST observations in the future can give more observational constraints on the gamma-ray 
emitting regions and the BLRs for the powerful blazars. 
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Fig. 1. — Sketch of axial cross section illustrating the geometry used for computation of the 
intensity of the BLR radiation and similar to Fig. 4 of Donea & Protheroe (2003). 
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Fig. 2. — The radial dependence of the energy density on the parameters tblr and h. The 
dotted curves in plots are calculated by adopting Lblr = 10 45 ' 35 erg s -1 , tblr = 0.43 pc, 
and tblr//cov — 1- From the bottom up, the four dotted curves have the parameter h at 
0.05, 0.15, 0.20, and 0.30 pc. The dash-dotted curves are calculated by adopting Lblr = 
10 45 ' 35 erg s -1 , h = 0.20 pc, and tblr// CO v = 1- From the top down, the four dash-dotted 
curves have the parameter t B lr at 0.30, 0.40, 0.50, and 0.60 pc. The plot (a) presents the 
energy density distributions for the emission lines, and the plot (b) presents the corresponding 
energy density distributions for the emission lines and continuum radiation. The plot (c) 
is the zoomed parts for the four dotted curves in the plot (b). The dotted vertical lines 
correspond to r blr = 0.43 pc, and the dash-dotted vertical lines correspond to tblr at 0.30, 
0.40, 0.50, and 0.60 pc. 
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Fig. 3. — The relative intensity of broad emission lines adopted by us. 
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Fig. 4. — The dependence of r 77 on the parameter h. In the calculations, we used tblr = 
0.43 pc and assumed i? 7 = r blr and t B lr// C ov = 1- The plot (a) is the photon-photon 
optical depth for four values of h, and the numbers attached to the four curves give h in 
pc. The plot (b) is the dependence of r 77 on h for four values of E 1 . The plot (c) is the 
normalized r 77 corresponding to the four cases in the plot (b). 
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Fig. 5. — The dependence of r 77 on the parameter i? 7 , and numbers attached to the curves 
give i? 7 in pc. In the calculations, tblr = 0.43 pc, h = 0.20 pc and tblr// CO v = 1 are 
assumed. 
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Fig. 6. — The dependence of r 77 on the parameter r blr> and numbers attached to the curves 
give tblr in pc. In the calculations, i? 7 = tblr, h = 0.20 pc, and tblr//cov = 1 are assumed. 



60 



50 



40 - 



30 - 



20 - 



l+ll 



10 



- 



50 



100 



150 



200 



E y (GeV) 



Fig. 7. — The optical depth of photon-photon absorption for FSRQs. The curve I corresponds 
to the absorption by the diluted blackbody radiation, the curve II the absorption by the 
emission lines, and the curve I+II the sum of the curves I and II. The three curves are 
calculated by adopting tblr = 0.43 pc, h = 0.20 pc, _R 7 = rBLR.in, and tblr//cov = 1- The 
curve III corresponds to the optical depth by adopting rBLR.in = 0.2 pc, rBLR,out = 0.6 pc, 
tblr/Zcov = 1/3, and i? 7 = r B LR,in- 
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Fig. 8.— The photon-photon optical depth for HFSRQ PKS 0405-123 are presented by the 
dotted curves, and numbers attached to the curves give h in pc. The dashed curve is the 
photon-photon optical depth for 3C 279, adopting r B LR,in = 0.1 pc and r BLR,out = 0.4 pc. In 
the calculations, i? 7 = r B LR,in an d t B lr/ f CO v = 1 are assumed. 



